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Abstract 


As the continuation of the previous report entitled “ Dynamo- 
theoretical conductivity and current in the ionosphere” (J.G.G. Vol. 4, 
No. 2, 1952), the effect of vertical motion of electrons on the electron 
density of F, layer is analytically treated. The vertical velocity of 
electrons in F; layer is connected to the density of horizontal dynamo 
current, which is obtained by solving the dynamo-theoretical equation, 
conserving the non-coincidence of the geographical and geomagnetic 
axes. 
The cyclic equation giving the electron density, wherein the 
term representing the effect of vertical motion is included, is solved 
with some approximations. 
The comparison of the calculated results with the observed data 
shows that the electron density depression and the height rise around 
noon near the magnetic equator, the longitudinal inequality of the 
latitude distribution of the noon density, and the shape of daytime 
variation of electron density for Ff, layer are plausibly interpreted by 
the present theory. 
I. Introduction 
In the previous report“ entitled “The effects of atmospheric motion and 
dynamo current on the electron density of the ionosphere” the author described that 
the horizontal atmospheric motion associated with the tidal circulation in the ionos- 
pheric level will affect the electron density of F, layer, but the effect of horizontal 
flow of charged particles corresponding to the dynamo current is much less than the 
former one. Then in the next report® entitled “ Dynamo-theoretical conductivity 
and current in the ionosphere ” the author investigated the ionospheric conductivity 
as of the heterogeneous and anisotropic medium, taking into account the so-called 
Hall-effect, and pointed out the significance of the vertical motion of the charged 
particles, which is produced accompanying the horizontal dynamo current. 
In the present report the author will describe some further extentions concern- 
ing the effect of the above-mentioned vertical flow of electrons on F, layer character- 
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istics, placing stress on the world-wide distribution of noon electron density and the 
daytime variation of F, electron density. The significance lies on the point that, 
while the former effect of the tidal circulation itself is based on the geographic coordi- 
nates, the present effect of vertical electron flow is associated with the geomagnetic 
coordinates and may provide a plausible interpretation to the geomagnetic distortion 
of the distribution and variation of F, electron density. 
II. Deduction of dynamo current on geomagnetic coodinates. 

As described in the second re- 
port®) we divide the ionosphere (F~ 
layer) into three zones. (see Fig. 1) 


O 


zone I; @=0°-75°, conductivity d,, 
zone II; 9@=75°-105°, 
conductivity 5d,, 
zone. III; @=105°-180°, ‘ 
conductivity Op. 
In this paper we use the same notation 
as in the previous one® and use@ and 


A as geomagnetic colatitude and long- 
Fig. 1—Division of zones I, II and III. itude respectively. 
Then the dynamo equations for the zones I and II are a follows, where the 
southward conductivity (dys) is considered as much greater than the eastward one 
(Oz) in the zone II, as discussed in the previous report. 
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and ¥,’, usually expressed with geographical coordinates, are expressed with geomagne- 
tic coordinates as follows.“ 


C5) Vy =R,'P1'(0) cos (A+ 44), 
=hy'{e1°P1° (9) +e,'P;! (0) sin(A+8;")}, 
(8) Yo? = ky? Po*(0) cos (24 +442), 


= hy? {€2°P2°(9) + cP 2'(0) sin (A+ Bo!) + €o2Po*(@) sin (2.4 + B2”)}, 
where 


¢,°=sin #)cos(T+A)+04') =0.2 cos(L+Ay +04"), 
2 (9) ¢i=V 1—sin? 6) cos? (T+ Ay +44!) = 0.98~1.0, 
£g81'=—cosf cot (T+A)+a,!)~—cot(T+Ay+ay5), 


€2° =3 sin*Oocos (2T + 249+ 002) = 0.12 cos (2T + 2A) + a2), 
- Ct=sind,V 1—sin?O, cos? (2T +2Ay+ a2) = 0.195 ~0.2, 


| — 
(10) a= 02 =y/costty}  Stcost (@T+ cos" eeu ae (2T +2 2+ 022) ~0.98, 
tgB.'!=—cosO)cot ears (2T +2Ao+ a2), 
tgp2= __1+c08"0o 04 (2T+ 2Ao +d?) ~—cot (27 +-2A9-F G2). 


2cosfy 


In the above, 4 and 4) are the geographical coordinates of the north geomagnetic 
pole, namely 6)=11.5° and 4,=—68°. We put ¢,', c.! and c,” as approximately const- — 
ant, namely 1.0, 0.2 and 1.0 respectively. Then the solutions of 7 or [, and J, are 
shown below. Putting 
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For the calculation of the right hand side of the above equation, I, and J, are re- 


placed by (12) and (13). Before the calculation 0 in the preceding section is described. 
Oy is the average value, which represents the Fy conductivity (common for southward 
and eastward) of the zone J, and 5d, is the conductivity (eastward) of the zone II. 
As for the vertical distribution of ¢;, we put oy as follows. 

e nN Vi 


m Ome | | Omi | 


(16) Op=Oyvz= CO = ie 


Taking the value of w,,, and w,,; as of 9=45°, and putting ADS MO, Ms, =O Tope 
the standard height, we obtain, 


n Ye Eg hs 
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Then (15) is calculated, using (17) and the equations of 6.2, Giz.....etc. and O’exy O'yy 


(17) Or=2 x 10° 


given in the previous paper.© 
1 0 f Vi fi | sin 9 cos ® nem me 
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(18) -- 2 (wa, = 


where fynZyn iS given in (12) and (138). 
In the above equations k,! and &,°, the coefficients of the velocity potential, are assumed 
as independent of height. And moreover the calculated results are valid for except near 


the geomagnetic poles. 


The first terms of 
the right hand sides of 
(18) are much less than 


the second terms, so far 


as the present considera- 


tions are limited to the 


low latitude zone or the 


daytime around local noon 


(cos@ and g’,, in (18) are 


very small.), and there- 


fore we can put appro- 


ximately as follows, neg- 
lecting the first term and 


ce 
OZNY-)- 


oS 
(19) (nv. ~—-K 5 


Geomagnetic colatitude @ in degrees 


Then K is approximately 
equal to v, and is given 
by the following equa- 


724 22 “20 “18 -16 “4 -/2 -10 -08 “06 ~04 -02 0 02 04 06 08 1.0 
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Fig. 2—Geomagnetic colatitude (@) versus ‘quantity proportional 
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3 sin® 
(20) Ke = ux Ries yn yn 


= 116% 10> ky nem gm, cm/sec. 


y was ae as 6,500km. In order to show the ieee latitude distribution of K, 
sin 0 DOEAS inn are illustrated for A=0° and 180° and for local noon (f=180°) in 
Fig. 2. In the calculation we put as hk,’ Fait , %1=0° and a.2=180°. K can be 
obtained by multiplying the value of sind 5314 a —*_ fm g™, (given in Fig. 2) with 1.16 
x10-°R',, that is, —3.48x10? (Rki!=—3.0x fpuasieh, In the present case, where we 
assumed (19), AK is identical to v,, that is, the average upward velocity of electrons 
in F, layer. From the curves we can see that the distribution of K is not symmet- 
rical for the north and south hemispheres (geomagnetic) and is also different for 
different longitudes (A). 
IV. The electron density of F, layer. 

The cyclic equation for the electron density (z) of F, layer is given as follows. 


(21) Ge q—ann—Bntrm Ko, 
where 7, and m, are the densities of the negative and positive ions. The strict solu- 
tion requires the set of three cyclic equations for m, m, and m, but as it is difficult 
to treat the problem strictly, we assume the terms containing a, f = ey hae aie? as 
follows. ~ 


(22) | 7=0n+B—-7p ( = 7), 


And also we Du “oF limiting our consideration within the daytime phenomena — 
near noon, which is assumed as steady. Then (21) becomes as follows. 


pan 
(23 cons pscll e 
(23) is Kae +7n=q. 
q and 7 are here given as the functions of z. AN prs 
(24) : ‘p= oe, > bs ; : get . : ‘G4 i as . * 


Distortional Characteristics of The World-Wide Distribution 89 


for K<0, 

m(0) pre fee e" 79 (e03 — gtr) 

1 u—'Sec x22" (ea eu 
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No 


- {(1—a)—sec x - 8 4s), 


The constant a in (24) is difficult to decide. @ must be less than unity, in order that 
the maximum of electron density should exist in the case of K=0. Considering the 
effects of an, and ys in (22) we assume the values of a and % roughly as 1/4 and 
3x 10-*sec™! respectively here in our numerical calculations. Then the equations (27) 
and (28) are numerically calculated for various values of K and x =0, putting H=50km. 

Fig. 3 shows the calculat- 


ed results of 2/nv,-»9 versus 3, 
where #-9 is the maximum 
electron density at K=0 and 
%=0, which appears at 3~—0.3. 


From these curves we can see 


the change of electron density 


associated with the vertical mo- 


tion of electron with the cons- 


oO or 02 03 o¢ Os 06 o7 og 09 40 
tant velocity. Electron density fatio (N/Nxag) 
Fig. 4 shows the two Fig. 3—Height distribution of electron density. 


groups of curves. The first group (solid line) shows the ratio of maximum electron 
density (z,,) for various K and x to m-o (the maximum at K=0 and x7=0). The 
second group (dotted line) shows 
the height (Zp) of z,, versus K 


Mee ok ey ee ; 
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(A) Geomagnetic latitude dis- 
tribution of noon electron density. 
According to the investiga- 

tion of Y. Aono® on the world-wide distribution of the noon electron density ratio 
(critical frequency ratio) of F, layer, which analysis was based on the world-wide 
data in 1946, 1947, 1948 and 1950, we can draw several relations of the noon density 
ratio (critical frequency ratio) as defined by Aono versus geomagnetic latitude 


+9 for various geomagnetic longitudes (A). Fig. 5 shows the three such rela- 


Fig. 4—Electron density and hp versus K (or v,). 
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Fig. 5—Relations between critical frequency ratio and geo- 
magnetic colatitude (plotted from Aono’s investigation). 
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tions for 1=0°, 180°, and 270° 
(90°) respectively. In the plot 
of the above relation we added 
several points (black points) in 
the case of 4=165°-195° of the 
bottom figure. These points 
are based on Japanese wartime 
(1941~1944) observations”? at 
fourteen stations from Kurile 
Islands (north) down to Bando- 
eng (south). As we are to dis- 


cuss here about equinoxial 
season, every point in Fig. 5 
represents the monthly median 
of one of the four months, that 
is, March, April, September and 
October. 

Against Fig. 5 the aes 


lated results are shown in Fig. 


6 of the ratio of the noon maxi- 


mum electron density (#,, at 
noon) to 2-9 (2, at K=0 and 
y=) 
obtained as follows. First given 
the longitude A, take several 
points of different colatitudes 
@ on this A, then the values of 
% at noon for these points are 
determined. Using Fig. 2, Kis 


These curves can be 


calculated for _ each point and | ) 
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Fig. 7—Daytime variation of critical frequency of 
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(Fig. 6) results. 

As the critical frequency is pro- 
portional to the square root of the 
electron density, one can see that the 
values near the equator of Figs. 5 and 
but for 


higher latitudes the coincidence is not 


6 are comparatively close, 


Satisfactory. 
(B) 
density. 


Daytime variation of electron 


Fig. 7 shows the daytime varia- 
tion of seasonal average (equinox) of 
Pv, observed at various stations. The 
data is partly based on Japanese war- 
time observations”? as described above, 
and the rest on observations at each 
country in 1947 and 1948. (A few data 
is observed in 1945, 1946 and 1950.) 
Most curves are the average of month- 
ly median of March, April, Septem- 
ber and October. 

Against these curves Fig. 8 shows 
the density variation for daytime calcu- 
lated for @ =90°, 75°, 70° and 60°, put- 
tinge "ay =25%). a2? = 204° Sand =A 180". 
As A and © are given, K is calculated 
for various times near noon from (20) 


72 
_ Local time 


43 


Fig. 8—Calculated daytime variation of - 
maximum electron density (F2 
layer) for A~180°. A curve of 
(cosy) ¥4 is given for comparison, 
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and using Fig. 4 we can draw Fig. 8. The agreement is pointed out below. 

(1) The existence of each one peak before and after noon near the geomagne- 
tic equator. But the afternoon peak at the equatorial zone locates a little earlier 
(around 15hr.) than the observed one (around 16hr.). 

(2) The occurrence of one peak shortly after noon at middle latitude. 

(C) Geomagnetic latitude distribution of layer height. 

Fig. 9 shows the average height of maximum ionization (Zp) at noon versus 


latitude (F-9). plotted from the observed results, which are mostly based on 
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Fig. 10—Calculated relation between height (Ap) 
of maximum electron density at noon 


for equinox and geomagnetic latitude 
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Fig. 9—Average height (kp) of maximum elect- Japanese wartime observations and 
ron density at noon for equinoxial partly on observations at each country 
pees versus geomagnetic latitude in 19pOeanal 195k 

—— ©) plotted f bserved data. ‘ , 5 
oe ) NaN eg Peach We can use Fig. 10, which is 


calculated for A~180° by using Figs. 2 and 4, for the comparison of our theory with 
the observed results (Fig. 9). The abrupt rise of height of about one hundred km 
for the low latitude is seen from Fig. 9 and the theory (Fig. 10) indicates the rise of 
about twice of 3, which is One hundred km with 50km of ZH. 

VI. Concluding remarks. 

The writer will discuss about a few points concerning the criticisms or ques- 
tions that may be raised against the validity of our present theory. 

(1) For such a solution of the problem as aimed at in this report it is a 
rigorous and orthodox approach to treat and solve simultaneously the following three 
equations. 

(i) j=f; (#, », 4) (dynamo equation), 

(ii) v= So, (% ¥; us J) (in short, for example, (15) of Section III), 

(iii) 2z=f, (q% v2) (cyclic equation), 
where the permanent earth magnetism and the velocity potential are given and 
assumed as independent of height, and y,; is assumed as given in a definite and con- 
stant proportion to v. When the height under consideration is fixed and q is given, © 
j, v, and m must be solved simultaneously from the above three equations. But such 
a treatment is very difficult to perform. The solution of (i) as a first step carried. | 
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out in our present study necessitated the fixation of conductivity, which can be defined 
by fixing 2, that must be the very goal to be reached at the final step. Thus our 
present theory has a contradiction in itself, although the writer presumes that there 
may be no substantial change in the results obtained. 

As the present theory indicates that the horizontal dynamo current produces 
the vertical flow of charged particles, and this vertical motion modifies the Chapman 
distribution of F, layer, it may presumably be possible to theoretically treat the 
ionospheric storm accompanying the magnetic disturbance. It will be the problem of 
future importance for the verification and extension of our theory to investigate the ° 
quantitative correspondence between the ionospheric disturbance and the disturbed 
current system associated with bay-type magnetic disturbance, or with the initial and 
main phases of magnetic storm. 
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On the Latitudinal Distribution of Maximum Values of 
f,F2 in Its Diurnal Variation 


By T. YONEZAWA and T. SHIMAZAKI 


(Radio Research Laboratories) 


Abstract 


Assuming that the temporal rates of electron production and 

removal in the F2 layer are of the forms q(cosz%)” and aN”, respec- 

tively, where x is the solar zenith distance, N is the electron density, 

and q, ¢, m and are constants, the latitudinal distribution of maxi- 

mum values of fyF2 in its diurnal variation has been considered, with 

the view of limiting the ranges of possible values for m and x. The 

results show that both m and mz should be contained in a range 

roughly between 1 and 1.5, which suggests that a process of attach- 

ment type must be predominant in diminishing electron density. 
1. Introduction : 
In considering the diurnal variation in electron density of the F2 layer, it is 
essential to know the functional forms which are taken by the temporal rates of 
electron production and removal. Usually, it has been assumed that the former is 


] 
} 
d 
z 


proportional to the cosine of the solar zenith distance x and the latter is of recom: | 


bination type, these being taken as qcosz and —aN”, respectively, where N is 


electron density, a is the effective recombination coefficient, and g is a constant. 


But the rate of electron production takes the above form only for an isothermal 
‘ 4 ¢ 


atmosphere at a level of maximum electron production. As it stands, the tempera: 
ture of the atmosphere changes with altitude and the height of maximum electr n 
production does not necessarily coincide with that of maximum electron £0 


ate = 


so. pes the rate of pelea gy generations na ee, a Pes a ges? tapers 
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limitation of the ranges of possible values for m and x», by considering the latitudinal 
distribution of maximum values of f)F2 in its diurnal variation. 
2. Method 

The method of analysis of data is as follows. The electron density variation 
in the daytime is governed by the following equation according to the assumption 


made in the previous section: 


Nis ewes ariouk> a) 
SF = qo (cos x)” —aN’ N-s lnT, (1) 


where ¢ is the time and T is the absolute temperature. On the right-hand side of 
this equation terms representing other effects such as the influence of winds in the 
upper atmosphere must be included, but assuming that these are of secondary 
importance,“ we have neglected them. At the time when the electron density 
becomes maximum in its diurnal variation, dN/dt=0 and, if it is not very far away 
from local noon, we may be allowed to take — lx T=0 also, so that we have at 
that time 
Jo(COS %»)" =AaN ms 


where x,, and N,, are the values taken by x and N, respectively, at the time of the 
maximum electron density. Now, considering that WN,, is proportional to f,,?, fin 
being the maximum value of f,F2 in its diurnal variation, and taking the logarithms 
of both sides of the above equation, we obtain 


log f,,= ae log COS X%»-+ const. (2) 


Strictly speaking, go and a in Eq. (1) will change with temperature T and T, in its 
turn, will depend on the solar zenith distance, and, consequently, the constant in Eq. 
(2) must be regarded as a function of cos7z,,. But, according to the estimation due 
to D.R. Bates,“ the upper atmospheric temperature does not seem to be much 
influenced by the solar zenith distance, so that we have here neglected this effect and 
assumed that the constant in Eq. (2) does not depend on z,,. 

Now, collecting as many as possible observed maximum values of f)F2 in its 
diurnal variation, f,,’s, and solar zenith distances at the time of maximum fpF'2, %»’s, 
at ionospheric stations widely distributed in the world, and plotting logf, against 
logcos x, we can obtain the value of m/2n from the slope of a line which fits to 
the distribution of the obtained group of points as closely as possible. In doing so, 
however, we must keep in mind that, while f,, may be dependent on factors other 
than cosz,, e.g. the geomagnetic latitude of the observing station or the height of the 
F2 layer, we must find out the variation in logf,, purely due to that in logcosz,, 
in order to obtain the value of m/2n, as is evident from the above exposition. For 
this purpose, we must subtract, at the outset, from logf,, the part which is due to 
other causes than 7, and the residual part must be considered anew as logf,, in 
: _ Eq. (2). In practice, taking the difference between the absolute values of geographical . 
; latitude, 2, and geomagnetic one, A, and the minimum virtual height of the F2 layer, 
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h’, as independent variables other than cos7z,,, and considering the multiple correla- 
tion between these variables and log/,, we can easily take out the part which may 
reasonably be deemed as purely depending on logcosz,,. On examining the data, 
however, h’ was found to have hardly correlation with f,, which allowed us to put 
it out of account from the beginning. 
3. Results 

Here we shall, as a first step, confine ourselves only to equinoctial seasons, and 
examine the data for March and September, 1947 to 1951, using Jonospheric Data 
published at C.R.P.L. The times of maximum f,F2 concentrate for the most part on 
12~13 hours, but in some cases they are shifted to times later than 15hour. Such 
an event must be brought about by some cause other than the equilibrium between 
the electron production by solar radiation and its removal by recombination and 
attachment, e.g. temperature effect or influence of winds, and as these effects are 
not taken into consideration in Eq. (2) and besides, as the points corresponding to 
these late maxima, when plotted, deviate markedly from the general tendency, we 
have excluded all these points.* The points corresponding to maxima before noon 
should have been excluded, when considered theoretically, for the same reason, but 


He as the arrangement of these 
16 


aly, we have adopted them. 


Aire sles alana A... RED, 2 SO a 
{I 


7 | = . appears near noon, in addition 
to the original one occurring not 


earlier than 16 hour, the former 
Sg ms has been taken. On looking at 

| ie ee 2 ‘ : A the scatter diagram of log f,, vs. 
Bs AG ks log cosx,, plotted according to 

the above-described criterion, it 
seems that there is no system- 


atic change in the distribution 


80 15 70° 
x 3 2 3 : 
, E met of points with the time of maxi- 
Fig. 1—Relation between the maximum value of f)F2 in ? 
its diurnal variation, f,, and the cosine of the mum, nor is there any marked 
solar zenith distance at-the time of maximum difference between the northern 


FoF 2, cosy%m, for March, 1951. Scales of both 
axes are logarithmic. Geomagnetic effect on 5 F 
foF2 is included. Pairs of letters attached to As an_ illustration we 


respective points are abbreviations for the loca- show in Fig.1 and Fig.2 the 
tions of stations. relation between logf, and 
log cos x,, for March, 1951. The abscissa is log cos %m and the ordinate is logf,,, but 
both axes are graduated in y,, and f,, respectively. The lines in these figures are 
determined to fit most closely to the scattered points in the least square sense. Fig. 1 


and the southern hemispheres. 


* In a few exceptional cases the maximum values of foF2 occurring at 15 hour are also omitted 
because they deviate from the general tendency too widely. 
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represents the case where f,, 
still contains the effect of geo- 
magnetic latitude and Fig. 2 the 
one where this effect has been 
removed from f,. One will 
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clearly see that the correlation 
rs “3 
between f,, and yx,, becomes BS} fy ADMaler) Mb 


i iG 
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more intimate in the latter s 
- . . — 
case. Figs.3 and 4 are similar ate Sa 
TY os 
figures for September, 1951. In 


this case the subtraction of geo- 


magnetic effect does not much 


improve the relation between 


: fo* . 60° ° yo" 30° 20700" 
Peeanay., for the partial cor: e hee Xe : es 
relation between /,, and | 2 | — Fig. 2—Relation between the maximum value of f,)F2 in 
| 4 | is low as shown in Table 1 its diurnal variation, f,, and the cosine of the 
below. solar zenith distance at the time of maximum 


FoF2, coSy¥m, for March, 1951. Scales of both 
In Table 1 are tabulated axes are logarithmic. Geomagnetic effect on 
the obtained values of multiple foF2 is excluded. Pairs of letters attached to 


correlation coefficient, partial respective points are abbreviations for the loca- 
> 


3 : tions of stations. 
correlation coefficient and x/m, 


and the number of ionospheric 


stations whose data have been 


used for calculation. Here f/f, x 
and A are abbreviations for 


log fm. log cos x,, and |A|—| A], 


respectively. 


Lastly, the most probable 


value of f,, can be expressed as 


a function of cosz,, and |A|—| A| 
in the following form which 
corresponds to the straight lines 
in Figs. 2 and 4: 

fo =FCOS Xn)" (LB)! 41-141 5 


(3) Fig. 3—Relation between the maximum value of JoF2 in 
its diurnal variation, fm, and the cosine of the 
solar zenith distance at the time of maximum 


here fo, a and 6 are constants 


whose values can be determined fof 2, cosym, for September, 1951. Scales of both 
i i i ic. Geomagnetic effect on foF2 
ideration of corre- axes are logarithmic 

peo is included. Pairs of letters attached to respec- 

lation between logfn» 108 COS Xm tive points are abbreviations for the locations of 

and |A|—|A| and are tabulated stations. 


in Table 2 for the respective months and years considered. In using these values in 
Eq. (3) .f, must be reckoned in Mc./s. and A and A in degrees. fo is the expected 
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Fig. 4—Relation between the maxi- 
mum value of f/)F2 in its 
diurnal variation, fm, and 
the cosine of the solar zen- 
ith distance at the time of 
maximum f)F2, cosym, for 
September, 1951. Scales of 
both axes are logarithmic. 
Geomagnetic effect on foF2 
is excluded. Pairs of letters 
attached to respective points 
are abbreviations for the 
locations of stations. 


Table 1. Values of Correlation Coefficients and »/m and Number of Stations. 


Month March September 
Year | 1947 1949 | 1950 | 1951 |Mean|| 1947 | 1948 | 1949 | 1950 | 1951 Mean d 
ryxa | 0.84 | 0.87 | 0.84 | 0.84 0.85 || 0.80 | 0.63 | 0.68 | 0.76 | 0.73 | 0.72 j 
rrxa | 0.77 | 0.78 | 0.77 | 0.82 0.79 | 0.75 | 0.44 | 0.68 | 0.74 | 0.71 | 0.66 
rrax | 0.55 | 0.59 | 0.63 0.60 || 0.57 | 0.49 
n/[m | 1.04 | 1.18 | 1.05 | 0.92 1.01 1.08 | 0.72 
Nun ber iatienliadinicl <cte|mersnienaeiatliletie Wars SO, al la 
of 41 41 
Stations 
Table 2. Values of Constants in Eq. (3) and Sunspot Number. : ; 
| Month || March | erneralek i aera 


|—________ 


| Year | 1947| 1948] 1949| 1950 | 1951 || 19 
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interest to examine these anom- isd 


16 


alies in detail, but here we & 
shall not discuss this problem 


further. Lat 


© Mar. 1949 


© Sep (947 
© Mar. 1950 © Sep. 1949 


From Table 1 we see that ° May. 1948 


4. Discussion of the Results 


the multiple correlation coef- 2 Sep 1348 


high (especially for March) so 


ficients, 7yx4’s, are considerably | 
that the obtained values of 2/m 


e@ Mar 1951 


may well be worthy of confi- ees 


dence. But it is a somewhat 

discouraging fact that the values o Sep 1951 

40 60 80 Joo 120 4o 160 180 
Sunspot Number 


of z/m cannot be considered as 
agreeing well with each other Fig. 5—Relation between the expected value of f,, for 
from year to year. On closer X¥m=0 and |%|=[A|, #, and sunspot number. 
examination we notice that this may be ascribed to different distributions of stations 
whose observational data are available for the respective years considered; each 
observatory has its peculiar characteristics and some are against general tendency, so 
that we may obtain unexpected results for years in which data from such perverse 
observatories are available in relatively large amount. Thus the dispersion of values 
of 2/m may be deemed to be inevitable to some extent at the present stage of 
observation, and we may be allowed to conclude, by adopting the average values, 
that 2/m takes the values of about 1.0 and 0.8.in March and September, respectively, 
without committing a great error. It is remarkable that z/m, in this way, takes 
relatively small values, for, according to the conventional method of treatment men- 
tioned at the beginning of § 1, m=1 and x=2, thus z/m becoming 2. 

If the upper atmosphere is not isothermal but increasing in temperature in the 
upward direction, 72 becomes larger than 1. According to M. Nicolet’s calculation, 
if the temperature (or, more exactly, the scale height) of the upper atmosphere is 


linearly increasing with height, 
m=1+f, (4) 


where # is the gradient of the scale height in the upward direction. Now, according 
to N.C. Gerson, it is possible that in the neighbourhood of the F2 region the tem- 
perature is increasing upwards by 500~1000° per 100km, so that, assuming that the 
atmosphere is composed of atomic oxygen, f, the gradient of the scale height, is found 
to be 0.26~0.52. Consequently, the upper limit to m may be considered as about 1.5 
by Eq. (4). On the other hand, the height of the level of maximum electron produc: 
tion rate expected theoretically may be considerably lower than the F2 region, being 
considered as roughly coincident with that of the Fl layer. If this is true, the optical 
depth of the level of maximum electron density in the F2 region will be much smaller 
than 1, and the rate of electron production at the level will not much depend on 
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the solar zenith distance, so that, if this rate is to be approximated by an expression 
of the form? q(cosz)”,m should be taken as considerably smaller than 1. Thus, even 
if the temperature gradient in the vicinity of the F2 layer is considerably steep, it 
seems possible that the value of m becomes much smaller than 1.5. As for the value 
of 2, since the main reactions conceivable as the mechanism of electron disappearance 
are recombination and attachment, we may be allowed to take m as not smaller than 
1 and not larger than 2. 

Combining these considerations concerning the values of m and z and the 


above-obtained result that 2/m=0.8~1.0, we can easily show that 
1.0<m<1.5, 1.0<r<1.5, (5) 


and that m and 2 are probably not very close to 1.5. Consequently, we might be 
justified to conclude that a process of attachment type is predominant in the mecha- 
nism of electron removal during the daytime also. But, as noted previously, the 
obtained values are not always well coincident with each other, and other seasons 
than the equinoxes are also to be studied, so that we should like to refrain from 
insisting on this conclusion with too much emphasis. 

Throughout the above considerations we have assumed the relation between 
logf,, and logcos x,, to be linear as expressed in Eq. (2), but this may not be the 
case. One who looks at Figs. 1-4 will feel that the distribution of points is not along 
a straight line but a curve which is somewhat concave upwards. We hope that we 
shall have an opportunity of discussing this point in a separate paper, and here we 
shall content ourselves with mentioning briefly that this fact implys that the slope, 
m/2n, in Eq. (2) becomes steeper and the constant in the same equation larger as %» 
approaches zero, and these, in their turn, might be regarded as evidences showing 
that the temperature gradient in the upper atmosphere becomes steeper as we come 
nearer the equator and local noon, and that @ in Eq. (1) decreases more rapidly with 
increase in temperature than does q in the same equation. 

In conclusion, some words will be spent on the partial correlation between 
logf,,, log cos x,, and |2|—| A 


. As is clear from Table 1, the high partial correlation 
between logf,, and logcos x, is a fact naturally to be expected, but the partial 
correlation between log f,, and |A|—| A| is, though consistently high in respective years 
for March, widely different from year to year for September, the difference being 
over the statistical fluctuation. In view of the large magnetic control over foF 2, this 
partial correlation is expected to be high, but the fact that this is not always the 
case may be the result of different distribution of stations whose data are available 
for respective months and years, like the case of the value of n/m, explained at the 
beginning of this section. In this respect it may be necessary to examine the charac- 
teristics of each observatory in detail, but this problem will not be considered here. 


5. Conclusions 


Assuming that the temporal rates of electron production and removal in the A 
layer are of the forms q(cos xz)" and aN’, respectively, the latitudinal distribution of 


—- 
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maximum values of f)F2 in its diurnal variation has been considered, with the view 
of limiting the ranges of possible values for m and nv. The results show that m and 
n should be contained in a range roughly between 1 and 1.5, which suggests that a 
process of attachment type must be predominant in diminishing electron density. 
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Self-Reversal of Thermo-Remanent Magnetism 
of Igneous Rocks (ID* 


By T. NAGATA, S. UYEDA, S. AKIMOTO and N. KAWAI 


(Geophysical Institute, Tokyo University) (Geological and Mineralogical 
Institute, Kyoto University) 


Abstract 


By separating the test sample of Haruna pumice at 300°C, the 
ferromagnetic mineral grains showing self-reversal of thermo-rema- 
nent magnetism (Reverse T.R.M.) (1) were classified into three types 
A, AB and B. The measurement of J,-T curve and T.R.M. character- 
istics of each type revealed that the sole component responsible for 
Reverse T.R.M. is the grains AB, each of which containing both of 
the constituents of types A and B in itself. It was concluded that 
the appearance of Reverse T.R.M. isa result of magnetic interaction 
between magnetic domains of the constituents A and B in the grains 


a 


—s 


of type AB during cooling under a weak magnetic field. 
1. Introduction pre “: 
In the previous paper (1), the writers reported an experimental proof of the 
self-reversal of thermo-remanent magnetism of igneous rocks, as well as of the 
separated ferromagnetic minerals from the rocks, and studied about the physical 
properties peculiar to this phenomenon. As to the physical mechanism of its causa- 
tion, two different possibilities which L. Néel (2) had proposed, were considered ; namely, 
the physical model producing Reverse T.R.M. as a result of magnetic interaction — 
between different constituents A and B, and the model containing N- or V-type | 
ferrimagnetic substance. In that discussion, it was also stated that, if the first possi- 
bility should be accepted, Néel’s original ‘model should be subject to some important 4 
_ modifications when applied to the present case. In other words, the writers s 
that the a pescnec) interaction must be within each individual grain and ‘the ¢ m1 
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the best way to separate out the constituent which is directly responsible for Reverse 
T.R.M. and to study various physical characteristics of that constituent. The experi- 
ments reported below are the products upon this view point. 

The experimental procedures consist of three parts: i.e. classification of the 
test sample by means of thermo-magnetic separation, measurement of the saturation 
magnetization and T.R.M. of each separated constituent. 

a) Thermo-magnetic separation of the test sample. 

Ordinarily, magnetic separation of minerals is done at room temperature. Since 
all the grains of the test sample are ferromagnetic at room temperature in the 
present case, magnetic separation at higher temperature was necessitated, which we 
call as thermo-magnetic separation. The method used was just simply to work a magnet 
on the specimen kept at desired temperature in nitrogen gas. Asa result, the separa- 
tion was inevitably limited to be only qualitative. Through magnetic separation at 300°C, 
it was proved that the ensemble of grains of the test sample which has previously 
been separated from the mother rock is classified into three types. By the reason 
explained in the next paragraph, we call these three types as types A, AB and B, 
respectively. The temperature of 300°C for thermo-magnetic separation was chosen, 
because the Curie point of the constituent B, T.2, was estimated to be below this 
temperature from the stepwise character of X-T curve (1) in a weak field and the 
Partial T.R.M. characteristics (1) of this sample. 

Although all of these three types seem undistinguishably ferromagnetic at room 
temperature, the circumstances are changed at 300°C, namely, though the grains of 
type A are still strongly ferromagnetic, those of type AB are feebler and those of 
type B are no more ferromagnetic. Although the quantitative results are still 
insufficient, the amount of the grains belonging to type B is much smaller than the 
others, as was expected from the fact that no traces of the existence of foreign 
constituents with lower Curie point were found in the J,-T curve of the sample as a 
whole (1). 

b) Saturation magnetization of separated constituents. 

As a next step, the J,-T characteristics of each type of grains were examined 
individually by a magnetic balance method. Taking into account the incompleteness 
of the thermo-magnetic separation, which will hinder the precise determination of the 
individual characteristics of each type, the measurement was performed on individual 
single grain. The results are shown in Fig. 1. As was expected, the J,-T curve of 
type A nearly coincides with that of ordinary magnetite or that of the test sample 


as a whole (1). On the other hand, type B shows extraordinary features in its J,-T 
en OJ, 
curve: i.e. the Curie point is about 230°C, | OT 


cerned temperature range, and the curve is linear, showing no tendency of usual 


is very large throughout the con- 


convexity even around the room temperature. From the result of chemical analysis 


of the test sample as a whole (1) and the extending work on magneto-chemistry of 


the system FeO-Fe,0;-Ti0: by E. Pouillard (3), this constituent may be considered to 


104 T. Nagata, S. Uyrepa, S. Aximoro and N. Kawar 


J be a kind of the substi- 
grkirary & i. tuted magnetites where 
some of Fet® of magne- 
tite are replaced by Tit# 
and Fe**. 

As to type AB, the 
forms of the J,-T curve 
below and above 230°C 
strikingly coincide with 
those of types B and A. 
Consequently, this type 


20,10; 5 


A | AB| B 
is regarded more likely 


as a mixture of the 


ret eet, aan T constituents A and B | 
Fig. 1—Changes of saturation magnetization of the grains of than as any other uniform | 
types A, AB and B. material. In other words, 
@ ....Heating process of type A 2 é 
n AB in- F 
© ....Cooling process of type A SE rate BS LE a ‘ 
— ....Heating process of type AB volves both of the cons- 7 
ia . Cooling process of type AB : tituents A and B within 
* eating process of type B itself. This isthe reason 
A. ...-.Cooling process of type B 


why this type was named 
as AB in the paragraph a). 

There was no indication of the existence of N- or V-type ferrimagnetic sub- 

; ‘stance in any of the separated types. 
c) Thermo-remanent magnetism of separated constituents. 

To verify and study in more detail the two constituents model, the T.R.M. 
characteristics of each type were examined. In this experiment, the measurement 
was also performed on individual single grain by which both effects of impurity. 
caused by the incompleteness of thermo-magnetic separation and the magnetic inter- _ 
_ actions among grains were perfectly avoided. Besides, the temperature range f field 3 
Bi ei re otocing, ILM. ae cise a by rok ee | 
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(3) Development of Reverse T.R.M. can be fulfilled by magnetic field cooling 


through the temperature range far below T,4, the domain fixing temperature of the 


constituent A. 


These facts obviously can not be explained by the original two constituents model. 


3. On the Mechanism of Reverse T.R.M. 


From the results obtained by the experiments so far, the physical mechanism 
of Reverse T.R.M. can be outlined so that all of the facts obtained hitherto are 
explained in a consistent way. 

The test sample is intrinsically composed of two constituents A and B of which 
the former will be nearly identical with ordinary magnetite Fe,O,; and the latter a 
kind of titanomagnetite, the solid solutions between magnetite and TiO.-2FeO = TiFe.O. 
If so, taking the fact that the Curie point of the constituent B is about 230°C into 
account, the content of TiFe,O,; in it is given to be 42% mol. according to the Pouil- 
lard’s result (3). The grains of types A and B consist of the constituents A and B 
alone respectively and those of type AB contain both of them. 

It is also ascertained that the grains belonging to type AB are the only res- 
ponsible component for Reverse T.R.M. and the mechanism of the modified two 
constituents model mentioned in the introduction of this paper should be realized 
within each grain of this type when the sample is cooled through the temperature 
range including the interval of R-process, say between 350°C and 200°C (1), under a 
weak magnetic field. According to this model, the various characteristics of this 
phenomenon, such as the Partial T.R.M. characteristic, the relationship between T.R.M. 
and the intensity of external magnetic field, and the development and: decay of T.R.M. 
can be consistently interpreted in almost the same manner as in the previous paper. 

The existence of type AB in the grains is to supply the more of evidence for 
the writers’ opinion stated in the previous paper that the magnetic interaction in 
question should be of domain scale. In other words, the said magnetic interaction 
would be in such the manner that the domains of constituents A and B are directly 
coupled so that the magnitude of magnetization effective for producing the demagne- 
tizing field would be of the order of spontaneous magnetization. 

With regard to the mechanism of the said magnetic interaction in the grains 
of type AB, there are three auxiliary experimental results supporting the writer’s 
inferences. The first is that Reverse T.R.M:s are produced through the field cooling 
in the temperature intervals of only 25° as long as they are included in the R-process 
range. The behaviours of development and decay with temperature were proved to 
be quite similar with the case of 50°-field cooling intervals. The second is that the 
domains of the constituent B in the type AB grains are fixed in the reverse direction 
even when the external magnetic field is strengthened up to 40 Oe. and the apparent 
T.RM. is in the normal direction. In other words, the effective demagnetizing field 
exserting on the constituent B is considered to be stronger than 40 Oe. in the case of 
the uppermost curve in Fig. 2, for example, because Jim 2 decreases with decrease in 
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temperature by the similar rate 
as in the case of Reverse T.R.M. 
The third is that the intensity 
of Reverse T.R.M. is decreased 
by grinding the test sample 
into finer dimensions. Jn the 
case when the mean diameter 
of the grains is diminished to 
less than 1 micron, it was ob- 
served that the intensity of 
Reverse T.R.M. J?),47-0.5 de- 
creases to some one tenth of 
the original case. This indicates 
that the structure of the grains 
of type AB may be destroyed 
mechanically. 

4. Conclutions and Further 
Problems 


x 
The mechanism of Re- j 
Fig. 2— 07 : i ot : 
‘ig. 2—Developments oF TR after cutting off the - verse T.R.M. in the present 1 
ternal magnetic field of various intensities. = 4 
The total process of field cooling is plotted for case seems to be clarified to = 
the case of H=0.50Oe. as an example. the extent of establishing the 
modified two constituents model. <= 


For further studies, the device for the more precise. thermo-magnetic cae 
tion is, now, under working. 

From the viewpoint of physics and mineralogy, the actual structure of the 
grains of type AB should be investigated by means of X-ray analysis, electron micro- 
Scope, mineralogical microscope or other. The nature of the constituent Bis also to 
be studied further. The measurement of its spontaneous magnetization, of which % 
change with temperature is extraordinary as shown in Fig. 1, to pee lower eee 
ture and other ER, Seabee reg may be of interest. ‘ Ries ~ . 


Bete aes > 
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York. It seems that more comprehensive studies are desired on these samples. 


(Read Oct. 6, 1952) 
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Development of S,-Field with Storm-Time 


By Takesi NAGATA and Hironori ONO 
(Geophysical Institute, Tokyo University) 


Abstract 


That the S,-field begins to grow up just after the sudden 
commencement of magnetic storm and it develops fairly intensely 
while the D,,-variation is still in its initial phase, is derived from a 
statistical examination of copies of H-, D-, and Z-traces of magneto- 
grams at Sitka (¢,,=60°.0). A simple graphical method to represent 
the variation of S,-field with storm-time is proposed. 
1. Introduction 
It was clearly pointed out by Chapman [1] and Vestine [2] that the magnetic 
storm field consists, in average, of the storm-time variation, D,:, and the disturbance 
daily variation, Sy. So far as Sp-field is defined to be the average deviation of daily 
variation on disturbed days from that on calm days, i.e. Sp=S,—S,, the distribution 
of S,-field may assume to depend only on latitude ¢, local time ¢ for different sea- 
sons. However, it has been well known at present [8, 4, 5] that the disturbance field 


% 
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quite similar to Sp-field begins to develop sooner or later after the beginning of 
magnetic storms. Then, the magnitude of S>,-field will also be a function of storm- 
time, T,. 

From the theoretical standpoint, on the other hand, the intense auroral zone 
current in the polar part of S»-field is attributed, now obviously, to the impinging of 
some charged corpuscular stream upon the upper atmosphere. Hence, the problem— 

“ At what stage of magnetic storm, the polar part of $3 field begins to take place ? Ei 
—may be an important key point for the corpuscular theory of magnetic storms. "7 
For the purpose of examining the above-mentioned problem, the copies of y 
magnetogram [6] at Sitka, Alaska (y,,=60°.0) during four years from 1946 to 1949 4 
were statistically analyzed, and the mode of development of S>-field with the storm- — 
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FeAl Sst); (1) 


where A(T,) denotes the amplitude of S,-field component at a point, while S(t) means 
the unit S,-field. If U.T. of the sudden commencement of magnetic storms and 
longitude of the station are denoted by T, and 4 respectively, there is the relation 
that 

=T,+(T,+4). (2) 


In the rectangular coordinate of ¢ in abscissa and T, in ordinate, (2) can be 
represented by a straight line making 2/4 with both coordinate axes and passing 
through ¢,=7T,+A on the abscissa. Then, the 
development of a magnetic disturbance with 
time can be expressed along the time scale 
on an oblique straight line in the ¢—T, plane. 
Different storms can be represented on different 
oblique lines. Consider now the distribution 


of disturbance forces F on the ¢—T, plane, SN 


on which a sufficient number of magnetic 


disturbances are plotted. Obviously, F is to 
be statistically shown as F(f,T,) on the f—T, 
plane. Further, provided that the assump- 
tion (1) is allowed, F(t,T,) is to be given 
by 


F(i,T,) = A(T;)- Sp). (3) 
The above-mentioned method of analy- 


Fig. 1 Current vector diagram of average 
sis was applied for H-, D- and Z-traces of the horizontal disturbing force of Sp- 


42 magnetic storms. It must be remarked component at Sitka. 


here that the polar magnetic storms, especially in and near the auroral zone, are com- 
posed generally of a number of elementary disturbances, the duration of each of 
which is from a few tens of minute to a few hours, as already pointed out by one 
of the writers and Fukushima [4]. Therefofe, F(é,Ts) concerned here are defined 
as those of the envelope-surface of a number of elementary disturbances. 

3. Results of analysis. 

First, the average Sp-variation at Sitka derived from all the used data is 
illustrated in Fig. 1 with vector arrow representation of the overhead electric current 
equivalent to the horizontal force disturbance. The distribution of the current arrows 
with respect to local time in the figure is well in agreement with the typical form 
of S,-variation [3, 4] and of bay-disturbance [7, 8] just outside the auroral zone. 
Since Sitka is situated at ¢,,=60°, the above result seems to be reasonable. The 
F(t,T,) diagram of H-component at Sitka is illustrated in Fig. 2, and those of H4D- 
and Z-components in Figs. 3 and 4 respectively. As obviously seen in these figures, 
the polar magnetic storm, which shows a sinusoidal variation in each component with 
respect to local time, begins to take place just after the commencement of storm, its 
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intensity increasing with the 


f 
| storm-time until 6~9 hours in 
ages pith 

Since S,(¢) is assumed to 
be a periodic function with 


SS 
~ H 
SS Se ar 


i Dee component of 4H at Honolulu 


12k respect to ¢, F(t,T,) in (8) can 
be written as 
é PV ROTI AB ATS sin| ty (t 
: / +4(T)}. (4) 
= a 6 12 3 18 ——-t ‘2h +50 0 —50 100" | The magnitudes of A; 


Fig. 2—Change in Sp-component in the horizontal (T,), As(T;) and #,°(T,) of AH 
* magnetic force with storm-time at Sitka. 
(£— Ts; Diagram) 


are plotted against 7, in Fig. 

5, where the average D,,-com- | 
ponent of 4H of the same | 
storms observed at Honolulu, 
Hawaii is also shown for com- 
parison. As will be seen in 
_this figure, both A,(7T,) and 
A;(T,) begins to grow up from 
0” in storm time, the former 
reaching the maximum value 
around 9%, and the latter 
about 5”. It will also be noticed 
that the phase (¢,°) of Sp-field’ 
with respect to ¢ is kept nearly 
constant throughout the whit 
‘Fig. 3—Change Cie ey inf AD onda atc 1~2h renee ; jst eine 
cig Ae Aye Te Deas et Say ig lina nee So 


te 
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zone starts to grow up from the SS + 


beginning of storm’s commencement, ot cote coin eek 
and it develops fairly intensely sale Ee rare eee 
while the D,,-variation is still in its Le 
initial phase. This fact may suggest an a 
that the corpuscular stream which is We ps 4 
responsible for the intense auroral zone stats Tis ae 
current of S,-field, arrives at the 1} [ates Seeeescriee aca 
upper atmosphere over the auroral i A2(Ts) | 
zone almost simultaneously with the 200|- ger ia San ae ae 
commencement of magnetic storms. a 

On the other hand, Fig. 4 being +54 im 
compared with Fig. 2, it may be hip clade’ (i ee aa 
generally said that the positive and - ee ia | 
negative values of 4Z correspond res- ar ae 
pectively to the positive and negative —_ 499 + 
values of 4H. This result will cor- [ l sais Tas e aaa 
respond to the wellknown fact that g : 2 ° i 15" 


—___» 


the main flow of the overhead auro- Fig. 5—Change with storm-time of amplitude 
ral zone current is generally situated and phase of diurnal component of 
at the north side of Sitka. Then, Sp-field in AH and that of amplitued 
the inversion of sign of 4Z about OSE eg Se ES eS 
7~10" in T, may indicate that the main auroral zone current is situated at the south- 
side of Sitka at the said time when the S,-field as well as D,,-field are in the most 
intense stage. This phenomenon may be due to the southward shifting of the northern 
auroral zone current due to the intensification of magnetic storm [9, 10, 11]. It may 
also be seen in Fig. 4 that the main auroral zone current returns toward north after 
the maximum activity of storms. 

4. Some Discussions. 

That the S,-field generally starts to grow up from the beginning of storm-time 
and becomes a complete form within about six hours in the low and middle latitudes 
has already been found by Chapman [12]. It seems that the present result concerning 
the S,-field near the auroral zone well harmonizes with the conclusion by Chapman. 
Hence, we may consider that the S>p-field over the globe begins almost simultaneously 
with commencement of storms and reaches its maximum intensity around the time 


of maximum of main phase of D,,-field. 
Another evidence for the above-mentioned conclusion will be the ‘simultaneity 


of the commencement of magnetic storm with the beginning of ionospheric storm in 
the F2-region in the auroral zone, which was found by one of the writers and Ogu- 
ti [13]. Since the ionospheric storms accompanying magnetic storms are believed to 
be caused by the impinging of corpuscular streams on the ionospheric region, the said 
result can be an evidence for the arrival of corpuscles at the ionosphere in the 
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auroral zone at the initial stage of magnetic storms. 

Let us now review theoretical pictures of physical mechanism of magnetic 
storms. It will be scarcely doubtful that the sudden commencement of magnetic 
storms is caused by the approaching of the front of corpuscular stream toward the 
earth, according to Chapman-Ferraro’s detailed investigations [14, 15]. Then, the 
fact that the auroral zone ionosphere begins to be disturbed soon after the time of 
sudden commencement at any local time may suggest that some corpuscles leak out 
from the main corpuscular stream which is retarded by the earth’s magnetic field, 
especially around the magnetic equatorial plane and those parts of corpuscles travel 
toward the auroral zone ionosphere along the lines of the earth’s magnetic force. In 
other words, it may be rather reasonable to assume that the front of the corpus- 
cular stream has already surrounded the earth at or soon after the time of sudden 
commencement. 

It will be worthwhile to note here that the S,-field in the initial phase of 
magnetic storm is not essentially different from that in the main phase and it is 
usually accompanied by the appearance of sporadic E-layer of about 100km in height 
in the auroral zone. This may mean that the energy of impinging corpuscles is 
about 0.3~1.0 MeV in its order of magnitude [16]. Whether the corpuscles have 
originally such the energy since they are emitted from the sun, or they are acce- 
lerated in the neighbourhood of the earth even in the early stage of magnetic storms ? 

This might also be a new question in physical theories of magnetic storms, 
together with the distribution characteristic of SC* with respect to latitude and local 


time [17]. These theoretical problems will be dealt with in detail in another article 
in the near future. 
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Influences of Ions in the Ionosphere* 
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Abstract 


It is examined that how the propagation of radio wave in the 
ionized region is influenced by the positive and negative ions in its 
region, and how a record is altered for the wave reflected from the 
ionized region for these ions. Especially, these probrems are discussed for 
the wave in and from the lower ionized region where the ionic density 
is considered to be great. The equation of motion and the Maxwell’s 
equation are solved in the presence of the earth’s magnetic field, ionic 
and electronic collisional frequencies, including the Lorentz polarization 

_ term. Results show that (1) influences by the two kinds of ions on the 
propagation become greater as the frequency of the wave decreases ; 
(2) energy of the reflection increases; (3) the deviation of energy and 
the virtual height of reflection for wave of certain frequency become 
smaller than those in the regions consisted of only electron; (4) the ‘ 
height of reflection thus approaches to the true height (no retardation), ; 
but the magnitude of the electron density, which is given by the 
relation N=A/?, appears to be doubtful in the case of the measure- 
. ment by low frequency radio wave. 
§ 1. Introduction ‘ : x 
The role which ions play in the ionized region has been olan by many 4 
workers in the study of the ionized region. This is because (1) it is considered that | | 
for the wave used in the observation of ionized region the wave is reflected from the | _ 
&€, F1 and F2 regions where the ionic density is small, therefore thew ave isa ffecte 3 
__ from ions and @) an ionic density is unknown in the stage of apr 7 J 


ee a ae! 
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small or great in the lower ionized region, hence we cannot easily determine. 

Recently observations on ionized region are often held by the waves of low and 
very low frequencies. We know that the heights where these waves are reflected 
are considerably low, though it is not known whether they are the tail of the 
FE region. The difference of the reflection height of these waves between day 
and night is abnormally large. It is, needless to say, admitted that the disapearance 
of electron at night is large for the great recombination coefficient at these levels. 
But besides, it is considered for us that the effect of ions may be introduced in the 
explanation for these results of measurements. 

In this paper we only discuss the role of ions on the propagation and reflection 
of wave in the ion-electron mixture. Its application to the phenomena in a real 
ionized region will be discussed in another paper. 

§ 2. Calculation of Complex Refractive Index 

The propagation of wave in a ionized region, consisted of positive and nega- 
tive ion and electron, was firstly discussed by Goubau.‘” He assumed that the positive 
and negative ion have all equal collisional frequency with other molecules and equal 
mass and density. Lorentz polarization term / was omitted without any explanation. 
In our study, we use same assumptions for ionic mass and density as Goubau did, but 
the Lorentz polarization term is not omitted in the equation of motion because we 
consider it is correct. This reason is explained in Appendix. 

We take z-axis in the vertical direction and x, y-axes in the east and north 
directions respectively. Vector of the earth’s magnetic field lies in the zy plane and 
makes angle a with z-axis. Direction of propagation of the wave lies along the z-axis 
and so the angle of incidence is zero. ; 


The equation of the motion for a iveed particle is given by 


P+ yp = Ne 


(1) 


where P, E, and H are vectors of polarization, electric field and earth’s magnetic field 
respectively; v is the collisional frequency of particle; N is the number density of 
particles per c.c.; m, e@ are the mas and charge of a particle; c is the velocity of light 
and, ‘represents the differentiation with the time. Maxwell’s equation is given for 


one kind of particle as follows; 


rot H=-(B+42¥Ne), | (2) 
rot B= —— H, (3) 


where r is the coodinate vector which has components 4%, y, z. For the mixture of 
electron and two sorts of ions 4zeNr is replaced by 4ndleNr. 
From (1). 

ee ‘ e ée 4 Gx e : as e 

x+ xv =a, E« + 4nlNex)+— (yH,—2H)), 


4 gy = © (Ey+-4nl Ney) + (— 8H), (4) 
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oo ° ée e 
24dy= © (B,-+4nlNez)+ 5 (8H). ) 
2, 
If E, x and y, z are proportional to e 5 bogae e ) then 
=“ [AE,+ Boosak,—sina-E,}, 
y =~ [—Beos aE, +(A—Csin’a)E,—Csina-cos E.], (5) 


2=—{B sing E,—Csina-cosaE,+(A—Ccos’a)E.]. 


with 
aS wv—v?—4n Ne?l/m ; 
(tov — w?— 47 Ne*l/m)*— ae o {iwv—w?—4n Ne?l/m} | tov —w*?—4r Ne?l/m)— i? E 
mc l mc 
nz 7weH/me ) 
(tov —w?—4r Ne hig a wo", Sa 
where w is the Feecoenee of wave, 4 is the wave-length, and yw is a refractive index. P 
q 


From (2) and (3) it is given in the ion-electron mixture 


HB = ee 2 ere = ies +4n31Nex], j 
H,= 9s ot = + (by +4n de SI, ; i ey 
| 1.=0, o= 1p. +4n5INed) , 
(7) leads to 
EAL —L)= —4rd)Nex, : 
E,1—p?)= —4n> Ney, hos Ig is (8). 
dey =—47)} Nez. : jes 
Substituting (6) into (7) and considering the gyro frequency ae shee eee 


ee shea) np pica ees in c com supe Sh ag Sattar: pesepoee es Ss “tag 
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ey 4nNe* cs 4nNie?* : : eH, a CH Ve 


Xe <a 9 wt yr = 
M,w* mw ? ~2" m.cw’ Iz Mm,CW’ A wo 


(suffixes e and z refer to ion and electron respectively) 


In the case of electron alone 2,2 is given by putting in (10) as #,2=1. 
§ 3. Reflection of Wave 

(a) Examples of dispersion curves in the propagation in the region of electron 
have been already given by Taylor. In the case of propagation in the ion-electron 
mixture, as can be seen from the expression (10), it is convenient firstly to draw the 
dispersion curve for propagation in the electron region and secondly to add the in- 
fluence of ions. It is clear that dispersion curves are similar to those in the electron 
alone. If Z, is not neglected the real of the refractive index does not come to zero. 
This tendency becomes greater as the collisional frequency increases for a certain 
wave. =a becomes zero when 


> 


Xe 


(l-)5=1-1Z,, (A) 

(i=) = SPR Syde (11) 
Rig a 

(I-)5=1-y-iZ,. ©) 


(A) corresponds to the ordinary wave and (B) and (C) correspond to the extra-ordinary 
waves. If “,;7=1, then these refer to the electron region. Further when Z,=0, (11) is 
expressed as follows: (Z is neglected) 


SS Pe 

m=15(1- 7H) AY 7 

Beis ae) ee (uy 
és rebige ep oes / 

%=15(1 : 4—\a eee ie 


If we consider that a reflection of the wave occurs when the refractive index is zero, 
even if x,=0, the reflection is possible when oe 

Generally, it is considered that the collisional frequency of electron cannot be 
neglected at the level where the ionic density is remarkably large (though the colli- 
sional freqnency of an ion is small in comparisom with that of electron and completely 
neglected), and so (11) cannot be, as it is, used at the level of E region or below the 
E region, but it is used in the Fl and F2 regions. For the magnitude of x; in latter 
case the concerning wave frequency is very large and x; does not easily exceeds 1 for 
the wave penetrated to these regions under the normal condition. For example, in 


Table 1 Table 1, orders of N; which makes x;=0.1 are shown. The wave 

y | Ni: that attains to the Fl and F2 regions is about 4mc/s (y=0.37) 
0.5 | 6.32108 and so we need 1.09x10°/c.c. as the order of N;. If the electron 
1.0 | 1.51108 density is 1.09x10%/c.c., 4 is about 10‘ and this is too great to ac- 


10 | 1.51106 cept. Even the wave reflected from the E region, N; should be 
100 | 1.51104 about 1.510% and x; does not exceeds 1 unless we adopt Booker- 
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Berkner’s aud Scott’s suggestion for 4. 

When y exceeds 10, order of Ni; which makes x; great abrupt decreases and 
therefore the effect of ions becomes greater for the lower frequency wave, which is 
reflected at the level far below the EF region. Consequently, we mainly need to study 
the effect of ions on the wave of low frequency. 

(b) When Z,-+0, real part of refractive index has not the value of zero or 
infinity. These examples are shown Fig. 1, 2, and for y=10, 100, Z,.=1, and 10 with 
imaginary parts without the influence of ions. Fig. 1 shows also curves including the 


acre WAVE 


EXTRAORDINARY WAVE 
y=10, 2.210 
ap=0 


Xe 
Fig. 1 Curves showing the real and imaginary Fig. 2 Curves showing the real and imaginary 
part of the complex refractive index, part of the complex refractive index, 
including the influence of ions for the including the influence of ions for the 
wave y=10. (Extraordinary wave). Hair wave y=10. (Ordinary wave) Hair lines 
lines and hair dotted lines show values and hair dotted lines show values of 
of R-(v) and [, (v) in ion-electron mix- R-(v) and In (y) respectively in the ion- 
ture. electron mixture. 
$00, 2210 ee influence of ions which has a uniform 


distribution through the region. Hair 
lines and dotted hair lines show real 
part of 4 (R,()) and imaginary part of 
“(Z,(2)). Full lines and dotted full 
lines show pens and ee mee 
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descends gradually to the point of the lower electron density than those in the case 
of electron alone. This implies that in the normal distribution of electron the height 
of reflection decreases. Further the absorption coefficient for a certain electron density 
increases because there the effect of ions adds. But considering the passage of the 
wave through the region 1t follows that the total absorption of the received wave 
will be remarkably diminished. As shown in Fig. 6, 7 and 8, the reflected energy in 


the propagation through the ionized region will increases. Variation of the reflection 


height of the wave when Z,-+0, is estimated from the following approximate expres- 
sion as 


3. 3—2%,  k=0 for ordinary wave 


= (lth) 53 +x; k=1 for extraordinary wave. oe) 


If the ionic region where the wave of low frequency is reflected has slow gradient for 
the electron density to the reflection level, the variation of reflection height reaches to 
the order of a few kilometer for etraordinary wave while for the ordinary wave this 
is neglected. For example, the extraordinary wave of y=10 is reflected at the level 
of 4500/c.c. of electron density and if this level lies at the maximum electron density 
and the region has the form of Chapman’s distribution which has the scale height of 
6km, then the variation of the height reaches about 9km in the case of fey hy [vene thi 
this region has maximum electron density of 10‘/c.c. the variation is about 2 km. 

Table 2 Hence for the wave reflected from the tail of the ionized region of 


Xe x; great electron density this variation is negligible. 


0.2 ; 0 ao. Table 2 shows the ratio of enough electron density of reflect 


0.4| 0.647 the incident wave in the presence of ions to the electron density in 
0.6 | 0.500 the absence of ions. As mentioned above, the influence of ions 
0.8 | 0.368 easily increases for the very long wave and so the electron deduced 


1.0 | 0.250 from the data of the observation of very long wave by the well- 


Pozi fences? known relation N=Af? will probably be incorrect and rather may 
4} 0.045 ete : : , 
: Bl 9 be greater that it is. Electron density assumed by Wilkes ef al., 


Bates and Seaton,“ Bremmer“ and Mitra“ must be reexamined. 


(c) As the ionic density is really 
not uniform through the ionized region, 0 
dispersion curves have a more compli- 
cated form. If we assume the ionic Fea(u 


Y=|0, Zél 
NASM 
4-0 


density is expressed as %,A AGA = Asean RARE 
fot. 2) and the electron density is 
distributed in a form of Chapman, then 


x; is given by 


x;= ene Nes eee) ee (13) Le 
eS Fig. 4 Variation of the real part of the refrac- 
For a certain 4», No, and a the disper- tive index for assumed o, # and electron 
12 


hi in Fig. 4. density distribution for the wave y=10, 
- sion curves are shown 1n a4 
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As can be seen in these figures, if x; varies sharply and has the maximum 
value, R,() dectreases abruply at the lower electron density than at the normal where 
the reflection occurs. In these points the abnormal reflection may occur, and if the 
considerable energy after a partial reflection can penetrate to the point of normal 
reflection, a second reflection occurs. If the first reflection is remarkably great, the 
second reflection will be unable to occur. If a partial reflection is possible, three or 
four reflection wave refered to the ordinary and extraordinary wave can be recorded. 

In order to examine whether the double 
reflection occur or not for one kind of 
wave it is convenient to use the *,—%; 
curves as shown in Fig. 5 in which 
dotted lines show the points of the 
normal reflection and full lines show 
the distribution of x; The reflection 
occurs at the neighbourhood of the 
_ intersection of two curves and even if 
- two curves do not intersect the reflec- 
tion occurs when the x; curve approa- 
ches near to the dotted one. There- 
fore if a curve firstly approaches the 
dotted curve and secondly intersects 
Fig. 5 The reflection point of the wave in the with this, then the double reflection is 
case of introduction of the ion. possible. To establish this condition 
it is necessary to examine the reflected energy. 
§ 4. Energy of Wave, Energy of Partial Reflection , ; 

The energy of wave is expressed by the real part part of the Seng Poynting 

vector hig) is given by 


E eas seis Hit (eH. nll, We-sFon 
Oe ae en ~ E,H Jk | 
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order to simplify (1 ; : d1Nez Nez 
plify (15) we make the following operations. As SiNex’? Sey are 
very complicated expressions and we must at first represent them in a simplified form. 


Before doing this it is convenient to start from the representation —_21Nez for the 


D>) Ne(x, y) 
electron alone. In this case 
d)Nez _ (2). >d\Nez Lee = 
aiNex (2). S\Ney (9),2 8 
where (x),, (y). and (z), represent components of r for electron. 
From (5) and (8) we get 
(Mf 1-11, 7 le 
16 
(shart (ZI tet 5)(9),=0, ee 
Lai pint Se =O 
From (16) we get 
: x 
: Z,—1—I x, £ 
(2)e —= Wr (2)e es — Me ; ee lene (17) 


(x). fc 4.—1—1x,+2,,° (We Wf, Vz 1Z,—1—lx, +x, ‘ 
Substituting (17) into (15), we get 


c zy 
(S.)e= — gy He(l— Me) |E.|? PVgES eee , 


Ce Aine Cer | ve 
(Sie = “go Hel — 7) y iZ.— 1- lx.— Xe tee 


Slee id Halt Fa) 


Further we get following relations 
Nee ee eens, Sen 
eae ree ere 2 oa 
rota I-23 yn yl Fett yet, yz} itd iz, y 


It is soe then that if we ner pe pollisrena) frequency ¢ of eon n the cel | 


(19). 


ae if ions exit inthe rion 
ERR IE 
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dWWez=(1+k)( Nez). 
hold approximately. Consequently it is given that 


DWWez _(@), DiNez _ ©), (21) 
DINex (x) Ney (y)." 


This implies teat influence of ions is approximately zero in the case of ratio. There- 


fore we can also use (17) in on ion-electron mixture. The influence of ions in 
the energy is introduced by replacing yu, by 4. Then the energy in the ion-electron 


mixture is given by real part of 


hee ea EE. 
S,=-—— u0-#)| Bala ta,” 
tz, eee fate 
fe : foe tte Lenk Ye (18) 
Sy=g, ed #')| Ey|? 1Z,.—1—Ix,—X- yr. 
Me poly 2 2 : 
S2= 3, ML | Ee | + | E,|?I. s 
If we neglect the collisional frequency we get 3 
: 
(Sede _Me(1— He") etc., (22) q 


Sek alae 


and if we introduce the collisional frequency we get 


R. (= x,2) Oil ’ 
aa d is ( a PY Eve oe 2She—x)dz (23) 
(Se Rie aie Zz oh ae. w Petey 
{( ZaLTe) 
where £,, is the initial amplitude of the incident wave in the x direction. The x, and x 
represents the absorption coefficient. The energy at any level in the ionized region is 
given by the addition of R,(S;), that is . 


RAS)=URAS)P HRS )PHRIS)PYE. (24) 


In this case, the energy reflected in the ionized region on the way till the primary 
wave attains to a pertain. eo is ace hence the variation of prince 


Saf 


we the Rabie distributed con a 
n e the 1 flection 
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where 2, 2, is the real part of “ in a narrow region 2 and 1, and Sha ©%1 is the 
(@) (aD) 


imaginary part of uw, %2, xX: is the absorption coefficient. The energy of the wave 
reflected from the upper boundary of the narrow region 1 is given by 


(72—1;)? 4 Gz | 
@ 


@ 
ors (26) 
(122+ 24)? + ( oo. oi) 


o wo. 


2R.(S):R=2R,(S) 


For simplicity, we assumed that that the energy reflected is not absorbed on the way 
to the receiver and also it is not reflected upward from the point below the level of 
the first reflection. The computation of R,.(S) is very labourious and so we take the 
special casees that there is no flow in x, y directions. These special places lie in the 
points of dip=0° and 90°. Then the energy is given by 


R(S)=R(S2)=G—-RAw) | Ee |2+| Ey) he? | (27) 


Y=100 , 20. 
—— Rel +) 


—| (S/Re( Sl 


1.0 


Fig. 6 Reflection coefficient of the incident Fig. 7™ Reflection coefficient, of the incident 
; 7 5 energy and reflected energy for ordinary 
energy and reflected energy for ordinary and extraordinary waves of y=100, when 


and extraordinary waves of y=10. Fos 


The ratio of the energy at a certain 
“et level to the energy of incident wave is 
can PA shown in Fig. 6, 7 and 9 and at the 
— R(S)/ RAS) ; : 
— RE same time the reflected energy is also 
shown. For a comparison, 4 for two 
kinds of waves are shown in the hair 


line. In these figures we assume that 


ae the electron density is distributed line- 
j x 100 HO ary and the thickness of the ionized 
Fig. 8 Reflection coefficient fof: the incident region is about 15km which is con- 
energy and reflected energy for ordinary sidered to be greater than the expected 
and extraordinary waves of y=100, when ; : ; 
Z-=10. thickness. From these domputations it 


is evident that the energy of the wave is reflected at a presumed electron density 
where the variation of the refractive index is large. For the ordinary wave in Fig. 
8 the explanation mentioned before does not hold and we must take another explana- 
tion that the absorption at a far distance from the normal reflection point, increases 
sharply as the wave penetrates in the ionized region and the energy reached there is 
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very small. It is too clear that the effect of ions is to increase the energy of the 
reflection. 

We can estimate by the same method the energy of the partial reflection for 
the then ordinary and extraordinary waves. The magnitude of the energy of reflec- 
tion depends on the point of the abrupt variation of the refractive index and the 
frequency of the wave. For example, if « decreases by 0.2 in the region of the 
breadth of semi-wave length as soon as the wave enter in the ionized region the 
sufficient energy to receive is given. This phenomenon may be possible in the practical 
observation and we may be able to receive the double reflection. But for the very 
long wave, as these reflection points approach each other it may be difficult to dis- 
tinguish the double reflection from the image of the receiver. 


§ 5. Deviation of Energy in the Ionized Region 
In connection with the flow of the energy in three directions, the deviation of 

the energy in the ionized region is an interesting problem. This is deduced from (18) 
or (18). If we take the ratio R(S,) to RAS.) and R.(S,) to R.(S,), the ratio of devia- 
tion of energy is given by 

R(S,) __—RAn(1— 1") Ey|?iyp/0Z,—lx, + 2— 1} 

RAS.) Rw Eel? + 1EyP | 

: A, 6 4 5 
BR (Sy) se) RAud—z’) |Eyl*y0(@-Z,— UXe+ X%_— | aE ) | ykid Ure eY) 


RAS.) F ROL | fe Pe | Be ae 


AeA neglected the deviation in the x direction is zero and (28) is changed to 


(28) 


RAS.) =2°)| Eg |) = bie) _ 22) iy, /{2(1-D 1} 


Re S;) (ORS 24 /2 2 
(S.) | a | Ey] tee tl 
: R. Sy = C= 1) yo( re Le— 1S Be \ [eet oe 
fe) ge 14142! ied = 
Further we take the following ‘special cases as dip=0° and 90, then a 
4 a5 ae = - ~ bh ES s mee . i eialk 
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the absence of the collisional frequency of electron. If the collisional frequency is not 
neglected this rule holds only soon after the incidence in the ionized region, but 
becomes invalid as the wave passes through the region and the deviation ratio changes 
the sign from + to — or from — to + depending on the magnitude of the electron 
collisional frequency, hence the energy flows often to the opposite direction. In Fig. 
9 the deviation ratio when Z, is zero is shown for y=1/2 and 10 in the region of 
electron. Curves for y=10 is different from those for y=1/2 and do not increase 
always as the electron density increases. Fig. 10 shows the deviation ratio in the 
presence of Z,. Two waves change the sign and especially the extraordinary wave 


— 4-05, 2-0 
---§=10 . 20) 4 


— E-WAVE 
--- 0-WAVE 


—s--—= ——— 


a 


Fig. 9 The magnitude of the deviation ratio Fig. 10 The magnitude of the deviation ratio 


and the direction of deviation of two and the direction of deviation of two 
components of the wave in the northern components of the wave in the north- — 
hemisphere in the absence of the colli- - ern hemisphere in the presence of the 
sional frequency of electron. collisional frequency of electron. 


twice. In the presence of ions this ratio for a certain electron density increases, but 
does not change the sign. Though the ratio of the deviation is large, the total devia- 
tion distance through the ionized region is rather small because the reflection height 
decreases. If the electron density is distributed linearly the area below the curve 
equals to the total deviation distance, and it is clear that the deviation becomes vx 
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where / is the coordinate of height. From (10) we obtain 


O/ 2) 372 
We wel 1— Sue} : (34) 
where ¢ is given by 
Abeer: gaa 
: eid “7 Own’ sin*a / Foto g'sin'a 
g=(o?+l0/u,?—toy,)— Se ara ot 3 7 Ree At ae 
w?— =) —twv, / (o- (1), —io»,) 
35 
If we differentiate 4 with w we obtain 2 
=( 6 \ a 
d E ne) dw 
Cer aie See 
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2 CHIEN ee ; : 
p= ees 1s shown in Fig. 11, 12 where it is know that the collisional frequency 


of electron influences rather than ions. Though p is greater a little for a certain x, 


In ion-electron mixture than that in electron alone, the difference of the virtual height 


is rather smaller. Especially, this tendency increases for large ionic density. Hence 


we regard that the virtual height of reflection of low frequency wave is in great 


measure correct. In the cases of dip=0° and 90° (37) represented by 


LES 2 )0-W2, 
aaa Reggc> Re ar ge eet | ordinary wave 


Maye 

° 4 9 bf 
yn | “3 4 «./ué—iZ,| 
(ISH Ore/ is — tale 


BiEXe 


(u."— w")2—1Z, 


=R| 1 + (39) 
extyaordinary wave 
dip=90° R, B32) R,| a+ u (2;°—w")(2—7Z, +1) for upper sign 
o A jiey. ordinary wave 
for lower sign extraordinary wave 


§ 7. Conclusion 
It is known that influences of ions become greater as the frequency of the wave 
decreases. By the existence of ions the absortion of wave decreases and the energy 
of the reflection become larger considerably than that in the case of electron alone. 
In the region of great number of ions, the estimation of the electron density must be 
reexamined. Deviation of the incident energy in the ionized region and the virtual 
height of reflection level is reduced. Especially, for a very long wave virtual height 


appraches closely to the true height. 
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APPENDIX 


On the Lorentz Polarization Term 


1. On the studies of Several Workers 

It is sutudied by several workers whether the Lorentz polarization term must be 
introduced or not in the equation of motion. Martyn and Munro(™ (M and M) 
assert the introduction of this term and Booker and Berkner“ (B and B) assert the 
introduction. Eckersley“ suggests that the phenomenon of existence of the critical 
frequency near 1.4mc/s is explained by the consideration that this occurs y;=1 and it 
is not important to persist the polarization term. Appleton, Farmer and Ratcliffe 
(A, F and R)“® also discussed this problem. Their conclusion are as follows ; 

(a) sort of wave causing a retardation | 

M and M=ordinary wave by a second reflection 

B and B=extraordinary wave 
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(b) region causing the retardation 

M and M=upper part of the F2 region 

B and B=below the £ region 

A, F and R=region between the EF, and FE, region or F2 region 

(c) magnitude of the frequency of wave 

B and B=Lorentz frequency 

M and M=f,cos0, where fy is the gyro-frequency and @ is the angle between 

the propagation of the wave and the vector of the earth’s magnentic 
field. 

We can point out the following defects in their studies ; 

(1) For the consideartion which B and B shows for the retarding region, M and 
show its defect. We also show the following point of view, that is, if the region con- 
sisted of a great number of ions, situates below the E region, the wave is reflected 
before it reaches the & region and a retardation does not occur. 

(2) M and M show that the wave is retarded at the level near the maximum 
electron density of the F2 region. But as the frequency of the wave incident in the 
F2 region is about 14mc/s it carinot at all reach such a height as M and M say, 
since the electron density necessary to reflect it is about 5x 10*/c.c.. . 

(3) M and M show the retardation occured by the ordinary wave, but for the 
ordinary wave the refractive index does not become infinity. Further they show that 


the frequency of retardation is f;,/cos0, but for the ordinary wave which they assert 
it deviation is northward in the northern hemisphere, then @ is about 90° and fycos#=0. 
(4) Eckersley says that retardation occurs at y,=1 from a fact that the retar- 
dation is not observed at lower latitnde. But the refractive index can become infinity 
near yp=1 and hence we must find another reason. ; 
_. (5) Though all workers miss the relation among y;, yp and x, (the electron 
density), it is important to consider their relation at different latitudes. 


— so — =-_ - 
J x 
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2. Mechanism of retardation Phenomenon 


In the explanation of the retardation phenomenon it is necessary to note five 
defects mentioned above and to know the behaviour of the extraordinary wave near 
(yz Or yp)=1. Some examples are shown in Fig. 13 and 14. 

It is easily found that the extraordinary wave behaves abnormally in a small 
range near (y,; or yr)=1. Collisional frepuency of electron also plays an important 
role there, and as the collisional frequency increases the maximum of R.() becomes 
smaller as can be seen in Fig. 13. 


If we neglect Z, the refractive index is given by 


jO=1 . é 
(SES ee SO ee eae Peas ye eee ee Tg eek & Rae aya’ 
es a=pel 92 =a) 4°" +y77(1—(1—2) x) f | 


For the infinity of the refracive index the expression is given by 


1 a 3 2 2 1 2 Zr 
(144%) (1-4.) —9e*(Lt+-g,)—915(1- Fx) =0. (41) 


Relation between y; and y, is easily determined by knowing the inclination and so 
from (41) the relation between y, or y, and x, (consequently y and x.) to make the 
refractive index infinity, is determined. This relation in the presence of / differs from 
that in the absence of J, that is, in the former the value of y exceeds generally 1 in 
the low or moderate electron density. On the other hand, in the latter, the value of 
y does never exceeds 1. In other words, in the fommer the wave frequency to cause 
the retardation is greater or smaller than the gyro-frequency in the retardaing region, 
and in the latter the wave frequency is greater than the gyro frequency. It is evident 
that the wave frequency of retardation . 


in two cases depends on the electron 


density of the retarding region. De- N eeoee \ 
Saulerare anon in dig. 15.) op ye ea 0 ‘ 
From the results mentioned  |0@@==s===-=~--__ tats. ' 


/ 
—_- 


above we explain the phenomenon of 


\ 
1 
Se XN ' 
critical frequency near 1.4mc/s as fol- fe * 
lows (we consider the case of dip=60°) ; ~ m 1 
(i) When the frequency of the wave se \ 
\ 


is lower than about 1.4mc/s, the ex- 


- 
t--- 


traordinary part of the wave does net 


- 
= 
ore 


cause the retardation in the EF region 


for the reason that the wave frequency 0 05 
is far below the gyro-frequency in the Fig. 15 Curves showing the relation between 
E region and it probably penetrates the y and x, which makes the refractive 
; index i i lies show the 
: tion of index infinity. Dotted 
Herecion though the aap , 5 case /=1/3 and full lines show the case 
this wave depends on the magnitude of ia: 


the electron density of the E region. 
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(ii) When the frequency of the wave is 1.38mc/s (Washington) or 1.4mc/s (Sydney), 
the wave reaches to the lower part of the F2 region (probably below 300km) with 
the electron density 5x10!/c.c. and there the wave is reflected in the absence of the 
special condition. But at this height the magnitude of the gyro-frequency is about 
1.43mc/s or more, hence from Fig. 20 the retardation occurs at. height with the elec- 
-tron density less than 5x10'/c.c. when we use /=1/3 (dip=60°) but does not occur 
when /=1. 

(iii) If the frequency of the wave increases to about 1.48mc/s the retardation of 
the wave in the E region also occurs. But there the absorption of the wave is very 
strong so that the wave is not clearly received. 

(iv) As mentioned above (ii) we are forced to accept the case 7=1/3 in order to 
accord with the observation. This leads to the result that the disappearance of the 
retardation in lower latitudes (about below dip=45°) is due to the fact that there the 
condition is similar to that of 7=0, that is, the retardation occurs when y<1. 

(v) As a result the frequency of the wave causing the retardathon is greater as 
the inclination increases. 

These retardation mechaism is different from that of other workers is the 
following point while simiar inpost exist. 

(i) The retardation is caused in the lower part of the F2 region under a certain 
condition derived from a relation between y and x,. This varies depending on the 
situation of the electron soa then the frequency of wave causing the retardation 
slightly varies. © 

(ii) The magnitude of the frequency of wave causing the retardation is different 
in each latitude for the reason that y varies considerably for different latitudes. 

(iii) The reason that the retardation of wave near 1.4mc/s is not observed is such 
that the value of y causing the retardation is smaller than unity in the lower latitude © 
below 45° of the inclination, and this corresponds to the case of 7=0. : 
(Read Oct. 5, 1952) 
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Propagation Characteristics of VHF over a Distance of 125km. 


By Kunio HIRAO, Yahei UESUGI and Kazuhiko TAO 


Abstract 


The statistical study of the field intensity of VHF wave over a 
distance of 125km was done by using mainly the observed data of 
150Mc wave and partly those of 65Mc wave. For this purpose, the 
weather maps, the aerological data obtained at Tateno Aerological 
Observatory and the meteorological elements observed at Kokubunji 
and Hiraiso were fully investigated. 

The field intensity and the range of its diurnal variation are 
large in summer and spring, and small in winter and on the day of 
precipitation. The abnormally high field intensity which appears 
frequently and its severe variation are considered to be reduced to 
the complicated distributions of the meteoroloical elements in the 
lower atmosphere. 

The effect of front is recongnized and it seems to be more 
remarkable in the night than in daytime. 

Introduction 

The statistical examination of the field intensity of ultra short wave has been 
done by many researchers, Report of the AD HOC Committee! is one of the most 
elaborate studies on the VHF propagation. But the period during which the data was 
obtained is not so long as covering a whole year. 

The routine observations of 65Mc and 150Mc waves between Kokubunji and 
Hiraiso have been continued since Dec. 1950. In this paper, these available data were 
used to make the statistical study for the ultra short wave propagation from the 
meteorological point of view. 

The Hiraiso Observatory is in 
the neinghbourhood of Mito and the 
we Kokubunji is at the westwards of To- 

kyo, and the path is across the plain 


"Fig. 1 Profil map of Kokubnnji-Hiraisa path. 
of Kwanto. The profil-map of this path is shown in Fig. 1. 

The present paper consists of the next five paragraphs ;— 
(I) Seasonal and diurnal variations of the field intensity and of its fading range. 
(Il) Aerological consideration of the field intensity. 
(II) Abnormally high field intensity. 
(IV) Effect of front on, the field intensity. 


(V) Conclusion. 
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I. Seasonal and diurnal variations of the field intensity and of its fading range. 


The fact that the diurnal variation of VHF wave is large at night and small in 
daytime, and the seasonal variation is large in antumm and small in winter, has 
been pointed out by many investigators since 1937.0 

Generally speaking, the above results seem to be correct, but there are many 
exception in each season on the days, when accompanied by special meteorological 
conditions. 

In this paper, the meteorlogical condition is classified into the following five 
types. 

i) Spring type; The travelling High covers the Japan Island. Warm and slight Wind. 

ii) Summer type; The Ogasawara High is situated at the Pacific southeast of Japan. 
Hot and slight wind. 

iii) Autumm type: Similar to that of spring. 

iv) Winter type: The Continental High covers the east side of Continent of Asia, 
and the ocean at east side of Japan is covered by cyclone. There 
is a large pressure gradient running from west to east over the 
Japan. Fine and cold, strong northwest wind. 

v) Type of the day of precipitation: The area of precipitation covers the whole 

— Japan Island. Rainy 
throughout a whole day. 

Those types have 
respective remarkable 
features for weather in 
Japan as ahown in Fig. 2. 

ee The mean diurnal 

' variations of the field 

intensity for every type 

: were calculated by aver- 
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Generally, the field intensity of VHF wave along this path is high in the night, 
especially in summer, low in daytime throughout all seasons. In the day of precipita- 
tion, the level is nearly the same all day long and in daytime nearly equal on all 
days. The diurnal variation in one group varies at the least in winter than in other 
seansons. This fact shows that the distribution of meteorological element in the lower 
atmosphere is nearly constant in winter. On the other hand, in summer it varies in 


considerable amount, showing that the lower atmosphere may be complicated. In the 


: next paragraph, it will 


be explained in detail. 
The statistical distri- 
butions of the level of 


[ 


a 


150Mc wave in each sea- 
son and that of a whole 


re 
c) 


year were calculated by 
considering the above 


a 


described meteorlogical 


& 


conditions. They are 
shown in Fig. 4. By exa- 
mining this figure the 
high field intensity seems 
to appear in spring and 


a. summer. In paragraph 
III, the anomalously high 
5 10 20 0 60 80 100 


field intensity will be 
discussed in detail. 
Further, the frequency distribution of the fading range of 150Mc wave at every 
five minutes were examined and it is concluded that the fading of a few decibels 
predominates in day-time, but in night  _ ["_ — 
time the fading range becomes large, 
especially during summer. The fre- 
quency distributions of fading ranges 
ecg ef aid winter types are 
? ee : 
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Fig. 4 Statistical Variation for 150Mc Wave. 
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Fig. 6 (A) The spectrum of fading range Fig. 6 (B) The spectrum of fading range 
of short period in summer. of short period in winter. 


examples of the vertical profil of refractive index are shown in Fig. 6. As known 
from this figure, there are many types of the vertical profil. But for brevity a linear 
distribution of refractive index will be assumed at first. The acutal occurence of 
the linear distribution for each month is tabulated in Table 2. This shows that the 
linear distribution of refractive index in the atmosphere appears more frequently in 
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Fig. 8 Frequency distribution of An, 
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cold season than in hot season. This fact coincides with that above described. 


Table 2 Frequency of the linear distribution of n. 
} Sy mains) or aod pe avi 
Jan. | |Feb. | Mar. | Apr. | May | June | July | Aug. | Sep. | Oct. | Nov. | Dec. rae 
Number of 31° 229g | <gpe {430 Mrste be 30-1 aa PS ABO eae 20 A sk eee 
observed days. 3 Sioa lio eae : ; Nt artes By an 
Number of Ui) 99.125. | oda) ba ahd 8: | 10S wk | Beil AO athe ue ieee 
% 94. 89 61 70 26 33 45 42 63 65 73 68 61 


The vertical gradient will be roughly expressed by the difference between the 
refractive index at the surface and that at the level of 500m high from the earth’s 
surface,—we denote this gradient by d4z—for convenience to compare the observing 
result with the theoretical one.“"1. The mean values of 4z for each month are repre- 
sented in Fig. 7. and the time of every value of dz is shown in Fig. 8. 4z in daytime 
are 16-20 MU and 20-30 MU in winter and summer respectively. The so-called equi- 


c— observed value 
wmumme Calculated value 


TFMAMITASOND 
Fig. 9 The monthly mean value of &. Fig. 10 The comparison of the observed field 
strength and the calculated one. 


valent radius of curvature of the earth, 
k, calculated from tha data in Fig. 8, is 


shown in Fig. 9 in which the level of 


k=4/3 is represented by a dotted line. 
The theoretically calculated value of 


the field strength estimated by dz in 
Fig. 7 is compared with the observed 
value in Fig. 10. The mode of propaga- 
tion during the day-time in in winter is 
so-called standard.‘ Within a whole 
year, the most frequent value of 4x is 
17 MU and that of the field intensity is 
5 db. 


As the theoretical field intensity 


Fig. 11 The relation between the meteorol- g P 
Giieat i eiuiiittons 2 uddl. the. Geli: is expected to be 4db calculated from 


strengths, _the value of 4m, the received intensity 
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is explained by the diffraction of wave around the earth which is surrounded by the 


medium that has linear vertical profil of refractive index. 

On the other hand, the theortical and observed field intensity were compared 
with each other for the individual days, two examples being shown in Fig. 11, where 
the theoretical value was not calculated, if the distribution of refractive index is not 
linear. Then, the observed values are generally high. This fact shows that the 
complicated meteorological condition causes the zig-zag distribution of refractive 
index and causes the anomalous propagation of ultra short wave. 


In Fig. 10, the simple descriptions of meteorological conditions are also entered 


as follows; d 
cf. ; cold front is in the neighbourhood of the path, 


wf 3; warm front is in the neighburhood of the path, 

H(a); travelling High has gone across the path, 

H(o); Ogasawara High spreads out near the path, 

L_; Low is in the neibourhood of the path, 

ws ; Low is over the Japan sea and the Plain of Kwanto is covered by the 
warm area. | 


III. Abnormally high field intensity. 


As mentioned in the previous paragraph, the abnormally high field intensity, 
which is unable to be explained by the linear distribution of the refractive index in 
the lower atmosphere, appears frequently. Although the meteorological process of 
the formation of the zig-zag distribution is not clear, the results in this paragraph 
will be useful in the future study. Fig. 12 shows a part of the period when the high 
field intensity more than 20 db of 150 Mc wave appeared during the period from Jan. 
1951 to May 1952. 

This figure clearly shows the results as follows; 

i) The abnormal high field intensity appeared more frequently in night time than 
in daytime. 
ii) When it once appears, it repeats in succeeding several days at about the cant 
time. ; 
In Fig. 13 the percentage of the 


ii Generally, the abnormal high 
intensity seems to appear when 
patente nee ise ca ie 


time of appearance of such a field to 
the total observed time is shown. This ; 
_ figure also shows the fact described in 
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intensity shows that, among the above six conditions, H(a) and H(o) is most effective 
to the anomalous propagation of- ultra short wave. Though the causual relation 
between them is not yet clear, the subsidnce may be possibly be the main cause of 
anomalous propagation in the cause of H(a) and H(o). 

IV. Effect of front on the field intensity. 

The propagation mode of ultra short wave, when the front is on the path of 
propagation, was investigatedlJ(9] and the observed effect was already pointed out.t71 

Now, the number of front which affects the propagation of the wave is examined. 
The effect of front on the propagation is remarkable in the field intensity and its 
variation. The former is abnormally high, and the latter is very large. 

The number of front which ‘passed over the path and the number of anomalous 
propagations which accompanied these fronts are tabulated in Table 3, and the per- 
centage of the latter to the former is shown in Fig. 14. The number and the per- 
centage are shown in Table 4 and in Fig. 15 respectively, where the cold front only 


is taken into consideration. 
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Fig. 14 The percentage of effective front. Fig. 15 The percentage of effective cold front. 


Table 3 Number of front and effective front. 


21-3 3-6 6-9 9-12 12-15 | 15-18 | 18-21 


Number! of front 36 22 22 25 24 45 33 


Number of effective front 20 9 9 8 3 la), 17 


Table 4 Number of cold front and effective cold front. 


Zio 3-6 6-9 6-12 12-15 | 15-18 | 18-21 


Number of cold front 12 12 Teh 13 is} 24 ae 


7 = 2 8 10 


Number of effective cold frot 15 if 


These figures show that the effect of front is little in daytime and that the 


effect of cold front is more predominant than that of the other type of front. What 


-kind of front affects the ultra short wave propagation is remained unsolved. To solve 


this problem, it is necessary to make clearer the structure of front, 
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V. Conclusion. 
It is concluded as follows: 
i) The field intensity of ultra short wave is high in warm and hot seasons, and is 


low in cold season and in the day of precipitation throughout a whole year. As 
for the diurnal variation it is considerablly regular, especially in winter season. 
ii) The structure of the lower atmosphere is analized, using the data of radio sonde. 
The vertical distribution of the refractive index is almost linear in winter but 
complicated in summer and spring. The field intensity on the day when the 
vertical distribution of refractive index is linear, is satisfactionally explained by 
the theoretical computation. 
iii) The abnormally high field intensity seems to ‘be in most cases caused by the 
travelling High and Ogasawara High. 
iv) The ultra short wave propagation is affected by front, especially by cold front, 
and the effect is more remarkable in the night than in daytime. 
In conclusion, the writers wish to express their thanks to Dr. H. Uyeda and to | 
Dr. T. Kono for their kind and instructive advice for the present study. They must 
pay also their thanks to the members of Tropospheric Propagation Section of Radio 
Research Laboratories, especially to the members of Unit of Radio Meteorology for 
their kind help in the work of the observation, calculation and others. 
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Ultra Short Wave Propagation and Mirage 


By Madoka FUKUSHIMA, Kazuhiko TAO, Yahei UYESUGI, 
Kiyoshi NISHIKORI and Kunio HIRAO 


Abstract 


The observations of the field intensity of 472 Mc, 153Mc and 
60 Mc waves and at the same time of the vertical distribution of the 
meteorological elements in the lower atmosphere were carried out to 
make claerer the anomalous propagation of VHF and UHF waves at 
the coast of Toyama Bay during a month in which mirage used to 
appear. 
On the day when mirage appeared, the field intensities of shorter 
waves increased considerably, especially for the 472 Mc wave. On the 
side of meteorology, the radio ducts were observed in these days and 
they were so intensitive as to transmit the 472 Mc wave in them. 
From this consideration, the mirage and the radio duct are thought 
to be reduced to the same meteorological conditions. The results of 
the present study will be available to forecasting the anomalous 
propagotion of ultra short wave in the neighbourhood of the coastal 
line. 
I. Introduction 
“Mirage” is a special phenomenon in meteorological optics. In Japan, it is 
frequently observed at coastal districts, especially along Toyama Bay during the months 
from May to July. It has been presumed that the air cooled by sea water is the 
main cause of “Mirage” at Toyama Bay, because the cold water of thaw of snow 
pours into the bay during this season and the temperature of sea surface becomes 
very low"), To affirm the fact that these meteorological conditions affect the ultra 
short wave propagation as same as the propagation of light, is the purpose of this 


experiment. 


II. Observation 
Both observations of the meteorological elements and of the field intensity of 


VHF and UHF waves were made at two towns of Uozu and Fushiki at coast of the 
Bay of Toyama. As to the observation of the meteorological elements, air tempera- 
ture and relative humidity, two types of theomometers and hygrometers, ordinary and 
electrical, were used. At sea shore at Uozu, a wooden tower, 12 meters high, was 
constructed and observation was made at the heights of 0.9m, 3.4m, and 113m from 


the mean sea level. At Fushiki a pole, 15m high, was used and observation was made 


at the heights of 0.8m, 3.4m and 13.7m from the mean sea level. 
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As to observation of ultra short wave, 60 Mc, 153 Mc and 472Mc waves, each 
horizontally polarized, were used. Among these three waves, 60 Mc wave was trans- 
mitted from Uozu and received at Fushiki and, on the contrary, other two waves 
were transmitted from Fushiki and received at Uozu. 

To magnify the effect of the anomalous propagation and made convenient to 

- compare the result with mirage, both transmitting and receiving antenna were set at 
lower level. Thus, the receiving point was situated at a point beyond horizon from 
the transmitting point in spite of short separation of about 29km. Conditions of j 


settings of these instruments are tabulated in Table 1. 


Table 1 Details of experiment. 


Frequency used 60.11 Mc 153.0 Mc 472Mc 
Transmitting antenna Height (above see level) 6.4m 15m 15m : 
ae ; (3.5m) : 
Receiving antenna Height (above see level) 20m 5.7m 3.4m 
(12.2m) ] 
Power 10Ww_ | ~ JoW 10W 


III. Results of observation. 

“ Mirages” were observed on 25th May and 2nd June. Time variations of 
modified index of refraction M at each level and the field intensities of electromagnetic 
waves are shown in Fig. 1 and 2. The time when mirages were observed at Uozu, is 
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Fig. 2 The observed data at 2nd June. 


shown in these figures by thick line. In every case, the duct appeared earlier at Uozu 
than at Fushiki and more remarkable at Uozu than at Fushiki. 

The field intensity of 472 Mc wave became anomalously high in these days and 
such a high field intensity didn’t occur on other days. Also, for 153 Mc wave, the 
intensity became considerably high in these two days at the same time. 

At to the 60Mc wave, the increase of the field intensity seems to be slight, 
although somewhat ambiguous due to interruption of the observation. 

At a glance, it is seen that the time of appearance of niirage, increase of field 
intensity and formation of a duct are not always coincident with each other. 

IV. Meteorological consideration. 

i) The weather map of the day of appearance of mirage is shown in Fig. 3. An 
anticyclone is on the sea eastwards of Japan and a part of it apreads out at the 
north part of Japan Sea and a local cyclone is at the central district of Japan. 
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Fig. 4 Mixing Ratio versus Hight Curve at 
Wajima Observatory. 


May 25th. 15.00, June 2nd 15.00 


Fig. 3 Weather maps of the days when 
mirage appeared. 


ii) The vertical distribution of mixing ratio of water vapour, calculated from the 
aerological data obtained at Wajima Meteorological Obsevatory, is repesented in 
Fig. 4. On May 25th and June 2nd, air at low level was dry; for example, the 
level with a mixing ratio of 5gm/kg is about 1600m and 1300m at 12h in both 
days respectively, while it is about 3000m at 12h in other days. Taking the 
pressure pattern in these days into consideration, the subsidence seems to have 
caused the dry air mass in the lower atmosphere. It was observed during the 
period of appearance of mirage. 
iii) At Uozu, it was almost fine and hot from early forenoon and a slight sea-breeze 
blew on the day mirage appeared. For several hours before the appearance of 
mirage, the air temperature became high and the humidity became low. Tem- 
perature increment and humidity decrement with height became so large that 
the difference of temperature and humidity between both observing points at top and 
the bottom reached about +3°C and —342, respectively, on June 2nd as extreme 
case. At that time, the temperature and relative humidity at the top were 28.8°C 
and 15% while the temperature of surface water was about 18°C. e 

In these meteorological conditions, a remarkable surface duct was formed and 
the ultra short wave trapped in the duct. At Fushiki, the weather was similar to 


that at Uozu. Ame: er —- wh ae nel ee 
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Fig. 5 Variation of mixing ratio at the day of mirage. 


inversion and humidity deficit. From meteorological point of view, temperature 
inversion and humidity deficit are both caused by the state of sea surface, so it 
is naturally considered that the mirage and the increase of field. intensity appear 
at the same time. The fact that is not the case, will be due to the non-uni- 
formity of the lower atmosphere. 

V. Consideration of the duct. 

It is obvious that a radio duct was formed in May 25th and June 2nd. The 
ducts are evaporation ducts those were formed by the above-described process. It is 
possible to estimate the thickness of these ducts by using meteorological data in a 
special case. In the present study, although the condition is not always satisfied, the 


approximate estimation was made by the following formula" ; 
7, 
d= P(q—9). 
3 


Here, d: thickness of a duct in ft, 
P: evaporation constant, 
qo: water vapour content at a standard level, 
q: water vapour content at a given level, 


72: 61x10! as 


fe 
73. 4./8— 59-px (14-46% 104.) per 100 ft. 
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The result of calculation is shown in 
Fig. 6, which shows that the height of 
the duct reaches 20 or 30 meters high 
during the most predominant period 
of the duct. On the other hand, the 
relation between the thinckness of the 
duct and the wave length of the wave 


which is able to trap in the duct, is 


Fig. 6 Variation of Duct Height. 


shown by the formula 
1=2.5dV 4M-10-* [4] 


The thickness of a duct in which 472 Mc wave is able to trap is about 20m in 
this case. The strong received intensity of 472 Mc wave accompanied by mirage, 
therefore, is explained as the field intensity of a superrefracted wave. As to the 153Mc 
wave, the field intensity may be more remarkable. 

VI. Conclusion. 

i) Mirage appeared on the day when meteorological elements had a special type 
with height. A hot and dry air mass travelled over a cold sea surface and a 
predominant temperature inversion occured in the lower atmosphere. This is 
the cause of the appearance of mirage. : 

ii) When mirage appeared, an extreme humidity deficit and a radio duct were : 
simultaneously formed in the same layer. This is the cause of the increase of 4 
the field intensity of the ultra short wave. Ew . 

In this experiment, the superrefraction of 427 Mc wave was observed, and the 1 
processes of the developments of these phenomena were studied. For other locations 
in Japan, where the similar meteorological condition appears, such anomalous propa- 
gation of VHF and UHF waves may be certainly forecasted. 

In conclusion, the writers wish to express their heartful thanks to Dr. AL Uyeda 
and Dr. T. Kono for their enouragement and instructive advice to the writers. 

This experiment and analysis were carried out by the members of the Tropos- 
pheric Propagation Section, Radio Reseach Laboratories. The writers wish also express — 
their thanks to them for their kind helps. eae 
“: The writers also pay their thanks to Mr. T. Ookubo and the mer 
a Police Station of apa 2 ease to the experim 
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